Spin-orbit splittings in heavy-light mesons and Dirac equation 



Riazuddir 

National Centre for Physics, 
Quaid-i-Azam University Campus, 
Islamabad 45320, Pakistan 

Sidra Shafiq 
Cenre for Advance Mathhematics and Physics, 
National University of Science and Technology, 
Islamabad, Pakistan 
(Dated: October 10, 2011) 

The spin-orbit splitting in heavy-light mesons is seen to be suppressed experimentally. It is 
shown that it can be understood qualitatively in the frame work of Dirac theory. An alternative 
derivation of a relativistic dynamical symmetry for the Dirac Hamiltonian, which suppresses spin 
orbit splitting, is also given. However it is shown that such a symmetry is not needed since the 
«4 ' spin-orbit splitting in Dirac theory with Coulomb like potential (as is the case for the one gluon 

I ■ exchange potential in pQCD) is small anyway. 

I. INTRODUCTION 

In the mass spectroscopy for heavy-light mesons as qQ or Qq bound states, one has both hyperfme splitting and 
spin-orbit splitting, which need to be understood in quantum chromodynamics (QCD). Writing a pseutoscalar (vector) 
heavy meson as P q (P*), P = D or B, q — u, d or s, the hyperfme splittings in MeV experimentally are [1] 

M D * - M D . = 140.64 ± 0.10 ~ M D , - M D * = 143.8 ± 0.04 

M B - ~M B - = 45.78 ± 0.35 ~ M B * - M B , = 46.5 ± 1.2 (1) 

"X' 

' On the other hand the spin orbit splittings seem to be suppressed (see below). 

C$ ' .... — 

We can write the effective Hamiltonian for a bound hadron containing one heavy quark(antiquark) Q(Q) as 

H = Hq + Hq (2) 

where Hq takes care of the residual momentum of the heavy quark (antiquark) and H q represent the motion of the 
light antiquark(quark) in a fixed potential provided by the heavy quark(antiquark). In heavy quark effective theory[2] 
(HQET), Hq contains a term q.~B* c /2itiq which gives rise to color magnetic moment interaction of type ~jl q .~jt q 
which induces the conventional form of the Fermi-Breit potential 

^a v (,){^^r) (3) 
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where ay(/i) = |ac(/x), ag is the pQCD running coupling constant which depends on the energy scale fx, in our case 

mass of the heavy quark. Here vhq and m q are effective constituent quark masses. The hyperfine splitting can be 

understood [3] in terms of the above term in Hq , which shows that such splittings decrease with the increasing mass 

of Q both because of toq in the denominator and decrease of a s (n) due to asymptotic freedom property of pQCD. 

— > — > 

Since in HQET the spin of the heavy quark is decoupled it is natural to combine j = L + S 9 , the angular 
momentum of light degrees of freedom, with S q to give J = j + S q for the bound Qq system. Thus one can have 
the following multiplets 

I = [P*(1-),P(0")] , 
I = 1 [P 3 *(2+),J\(1+)] 3 

J J— 2 

[P 1 *(1+),P (0 + )] J=| (4) 

where P is D or B and J p gives the total angular momentum and parity quantum numbers. The splitting between 
3 = § an d j = \ for I = 1 is due to the spin-orbit coupling L .S q while the hyperfine splitting between 2 members of 
each multiplet arise from the Fermi-Breit term as mentioned earlier. The spin-orbit splitting (in MeV)[l] for D mesons, 
between D* 2 {2+) : 2462.8 ± 1.0 and : 2422.3 ± 0.6MeV is 40, for the B mesons, between P^(2+) : 5743.9 ±5.0 

and Bl : 5723.4 ± 2.0 is 21, for Bf 2 {2+) : 5839.7 ± 0.6 and P Sl (l+) : 5829.4 ± 0.7 is 11. Thus these splittings are 
suppressed. A measure of this suppression is the parameter [4] 

P3/2-P1/2 



(4p 3 /2 + 2pi/ 2 )/6 - 

which for the experimental data shown above is of order 0.07 both for D and B mesons. 



(5) 



II. A DYNAMICAL SPIN AND ORBITAL ANGULAR MOMENTUM SYMMETRY FOR THE DIRAC 

HAMILTONIAN 



As it is well known, the solution of the Dirac equation with Coulomb potential for the hydrogen atom does give 
spin-orbit splitting between 2ps and 2pi energy levels in agreement with experiment. Since the one gluon exchange 
potential (OGE) in pQCD is Coulomb like, one would expect [5] such a splitting in the spectrum of hadrons, in 
particular heavy mesons considered here. Thus we take H q in the Eq. ([2]) as the Dirac Hamiltonian [setting h = c = 1] 

H = M Q + ^.-f + p{m + V s ) + V v (6) 

where ~^ — —i^ is the 3-momentum operator, and j3 are the Dirac matrices, m is the mass of the light quark 
and Mq that of the heavy quark Q. The above Hamiltonian as remarked earlier describes the motion of light 
quark(antiquark) in a fixed potential produced by the heavy antiquark(quark) as in the hydrogen atom where electron 
moves in Coulomb potential provided by the nucleus (proton). We have assumed that vector and scalar potentials are 
present, the latter for the reason to be stated shortly. 

It has been observed [4,6] that if vector and scalar potentials satisfy the relation 

V v (?) = V a (l>) + U (7) 
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where U is independent of the position of the light quark relative to the heavy quark, then the Dirac Hamiltonian is 
invariant under a spin symmetry (called relativistic spin symmetry) 

[H,Si]=0. (8) 

and if potentials are spherically symmetric, then there is an additional symmetry 

[H,Li]=0 (9) 

The generators of these symmetries are given by 




here Sj = ^ are usual spin generators, <7j the Pauli matrices, ij = (r x p)j while Sj = U p SiU p and l, L = U p liU p , with 
U p = ^jp as the helicity operator. Thus even though the system may be highly relativistic, the Dirac eigenstates can 
be labeled with orbital angular momentum as well as spin, and the states with the same orbital angular momentum 
are degenerate, e.g the states n r pi and n r p3 are degenerate where n r is the radial quantum number [4]. 

Thus a symmetry has been identified in the heavy-light quark system which produces spin-orbit degeneracies 
independent of the details of the potential. It may be pointed out that for the hydrogen atom, such a symmetry 
would be contrary to experiment since the splitting between energy levels 2p 3 / 2 and 2p 1 / 2 as predicted by Dirac 
equation is in agreement with experiment which has been regarded as a great success of Dirac equation. The levels 
S1/2 and pi/2 are still degenerate (the Lamb shift) which requires quantum radiative corrections. 



III. DIRAC EQUATION AND SPIN SYMMETRY 



By writing the Dirac equation in two component Pauli form, the symmetry discussed in the previous section can 
be derived in a much more transparent way which would also help us to solve exactly the Dirac equation if both the 
vector and scalar potentials are Coulomb like. Consider the Dirac equation in covariant form in the presence of a 
gauge field 

(i^Df, - m)tf = (10) 

where = 8^ + iV^ and 7 = (3, 7' = fid 1 . For our case taking the zeroth component of and a scalar potential 
V s , the above equation become 

[i^ido + iV v (r)) + i 1 l d l to V a (r)]tf = 0. (11) 

Multiplying on the left by 

[if{d + iV v {r)) + i^dj + m + V.(r)] (12) 
we obtain, since the potentials are independent of time, 

+ 2iV v (r)d + V v 2 (r) - + nVt^, V v ] ~ if{d u V s ] - (to + V s (r)) 2 }* = 0. 
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For stationary states, % — > iE, and using the fact that djdi is symmetric in j and i, we obtain 

[V 2 + V 2 - V 2 - 2EV V ~ 2mV s + H°f[d t , V v ] ~ if^, V s ] + (E 2 - m 2 )]* = 0, (13) 



where 



and the second equality holds for spherically symmetric potentials. 
It is convenient to use the chiral representation of 7— matrices, 



7° = P = 



>7 5 = nV7 2 7 3 = 
10/ \ 1 



, -o-' , 

a* = /3 7 ' = ■ (15) 



cr' 

Then we can write the above equation in two component matrix form 



01/ \ / / \ 



where 

O = V 2 + V 2 - V 2 - 2EV V - 2mV s + (E 2 - m 2 ). (17) 

The diagonalization of the matrix multiplying icr.f , since (cr.f) 2 = 1, gives the eigenvalues rj 

dV v dV v ,dV s 2 n 

or or ar 



?? = ±[( ^ )2 _(^ )2 ]i (i 8) 

The matrix which diagonalizes this matrix does not affect the first term in the equation (|16[) as it is multiplied by a 
unit matrix which commutes with every matrix. Denoting the corresponding eigenfunctions by ^ ± which are linear 
combinations of ^ l and ^n, we have 



[6 ± ir)(T.r}V± = 0. (19) 



We note that the eigenvalues 77 in Eq. (fT5|) vanish for 



dV v dV s , x 

IF = ± (20) 



V v (r) = ±V s (r) + constant. (21) 
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Then 

6tf± = 0, (22) 

where O is independent of spin. If plus sign is selected, the system is said to have relativistic spin-orbital angular 
momentum symmetry; for negative sign it is known as the pseudo spin symmetry[6] which has been observed in 
nuclei [7]. As a result there is no spin-orbit coupling and the results obtained in the previous section are derived in a 
different and more transparent way. Selecting the positive sign in Eqs. (|2ip and (|22p we can solve the equation 

[6 ± iria.r]$± = 0, (23) 

exactly for the energy eigenvalues as in hydrogen atom for 

V v = -^ + U V , 
r 

V s = -— + U S . (24) 
r 

Here V v (r) is the vector OGE potential with a v — |q;g in which we are interested and V s (r) is the confining scalar 
potential which might arise from multi gluon effect. But its origin to be Coulomb like is not clear (see, however ref 
[4] and references there in). 
For the above case we have 

(a 2 -a 2 )5 

r) = [ s ^ . (25) 

Using the spherical polar coordinates the operator O becomes 

A d 2 2d L 2 a 2 — a 2 n/ i > .1 . „. 

° = 7^ + Z-- — + + 2 (£ + m a s )- + (E 2 -m 2 ), (26) 



where 



Accordingly the equation (|19l) becomes 



E = E-U v , (27) 
m = m + U s . 



(a 2 -a 2 )^ 

[0 + i { v " s) a.T]^ + =0. (28) 



The energy eigenstates can be now read of from those for the hydrogen atom [8] by making the substitutions 

E' , , 

a — > a v — 7 +a s , (29) 



rn 

E 

E -> 



E 2 - m" 2 



2m' 

and are given by 



2m' 2{n-5,j) 2 y ' 
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FIG. 1: Matching of OGE (Coulomb like potential, dashed) with Cornell potential (Coulomb plus linear r 

dependence, solid) 



Solving this quadratic equation for E' and re-substituting the values of primed quantities given in Eq. (1271) 

E n j = U V + -ZTWZ 2 ["Q^s ±(n - 5 )[{n - Sj) 2 + a 2 v - a 2 s ]^}, (31) 

where 5j is 

Sj=(j±l)-{(j±l) 2 -(a 2 v -a 2 s )^. (32) 

This is an exact result and n is the principle quantum number. When V v — V s , so that a v — a s and U v = U s 5j 
becomes zero and 

2a 2 

As the expression for energy is independent of j, there is no spin orbit-splitting for this case. 



IV. SPIN-ORBIT SPLITTING AND CONCLUSIONS 



We now discuss spin-orbit splitting based on equations (l3Tj) and (|32|) . For our case 1 = 1, n — 2 and we have 1 + 
and 2 + mesons corresponding to j = h and j = |. Since experimentally 2 + mesons are heavier than 1 + mesons, i.e. 
E 2 3 > E 2 i , one may conclude that Si > S3. Then Eq. ([3^)) implies that (a 2 — a 2 ) > and further since Si is real 
(a 2 - a 2 ) < 1. For V v {r) we take OGE potential 

K(r) = -— + U V (33) 
r 

where a v is treated as an effective coupling constant. We put a s — and represent the scalar potential by a constant 
U s , the whole purpose of which is to renormalize the light quark mass in hadron to m + U s [c.f Eq. pip ] where U s is 
to be fixed by the data on spin-orbit splitting. How good is this approach? To see this, for the charmed sector, we 
compare the potential given in Eq. (j33| with the Cornell potential [9] 

V{r) = -— + ^ + C (34) 
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with K — 0.48, a = 2.34GeV r ~ 1 and C — —Q.25GeV. This comparison is shown in Fig 1 with the boundary condition 
(/ stands for fermi), T^(0.5/) = and the matching point r = 0.14/, which gives a v = 0.8, U v = 0.32GeV. With 
these two conditions, which the Cornell potential also satisfies, it is matched with the lattice QCD potential[10] in ref. 
[11], showing that the Cornell potential gives the simplest extrapolation of the lattice QCD potential. The potential 
given in Eq. ([33)) with a v — 0.8 and U r = 0.32GeV almost give the same extrapolation as is clear from Fig.l. 

We now proceed with the numerical results. First we note that the mass of Qq or qQ meson is given by (n = 2) 
[E 2 j is given in Eq. (j3"Tj) ] 

M 2j =M q + E 2j . (35) 

Then the mean mass of j = A, j = | mesons is 

M = M Q + i(£ 2 , f + E 2>h ) = M Q + U v + \{m + U.)(F 2 , + F 2 ±) (36) 
while the mass splitting is given by 

AAf - (£ 2if - E 2 i) - (m + U S )(F 2 ^ 2 - F 2>| ) (37) 

where [with a s = 0] 

^ = ro_ A + „2 [( 2 - J i)[( 2 - + a ^ ( 38 ) 

To carry out the numerical work we have to fix light quark u or s mass m, Mq and U s . It is known from the mass 
spectra of I = mesons that 

m = m u = 330MeV 

m s = 55QMeV 

M c = MSOiWW 

M b = 4800^/6^ 
4 

a v = g"G, (39) 

where a v and U v have already been fixed for the charmed sector. As the QCD coupling a v is energy dependent, 
decreasing with increasing energy, one would expect 

a v {M B J < a v (M B ) < a v (M D ), 

we will take this into consideration. 

It is clear from Eqs. (|3"6")l , (f3"Tf that once we have calculated F 2 j, (m + U s ) can be fixed from Eq. (|3l))) . Our numerical 
results are summarized in Tables I and II. 

To conclude we see that the spin-orbit splittings, obtained from the Dirac equation with OGE potential and scalar 
confining potential to be a constant so that its role is to enhance the mass of light quark in a hadron, are small and 
qualitatively explain the data within a factor of about 2 (for charm) and 1.5 (for bottom). The values used for a v 
for the OGE potential are some what larger from those obtained from the asymptotic freedom of QCD. However, in 
potential models a v is usually treated as a phenomenological parameter. From Table 2, we see that U s for the scalar 
confining potential is almost independent of flavor. Another important conclusion one can draw is that as in hydrogen 
atom, the meson spectroscopy does not show relativistic spin and orbital angular momentum symmetry, since the 
splittings are small anyway. 
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TABLE I: Experimental data; all masses in MeV 







Mesons 


Mass M 


Mean Mass M 


Mass Splitthij 


y AM 










2422 


2442 


40 










D* 2 {2+) 


2462 














Bt{\+) 


5723 


5733 


21 










Bi{2+) 


5744 














B* S1 {1 + ) 


5829 


5834 


11 










b; 2 (2+) 


5840 
















TABLE II: Predicted mass splitting 


;s m MeV 






Flavor 


a v 


u v 


F 3 - Fi 

■?. 9. 


(F 3 + Fi)/2 


m + Us 


AM(our) 


AM(expt) 


D 


0.8 


320 


0.029 


0.901 


712 


21 


40 


B 


0.75 


300 


0.019 


0.917 


690 


13 


21 


B s 


0.65 


260 


0.009 


0.938 


814 


7.3 
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